Journal of the College of Physicians and Surgeons Pakistan
ISSN: 1022-386X (PRINT)
ISSN: 1681-7168 (ONLINE)
Affiliations
doi: 10.29271/jcpsp.2024.11.1831ABSTRACT
Objective: To construct and validate a prognostic model for osteosarcoma prognostication and therapeutic potential of PANoptosis- related genes.
Study Design: Observational study.
Place and Duration of the Study: Department of Orthopaedics, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China, from August 2021 to January 2024.
Methodology: Transcriptomic data from the GEO and TARGET databases were utilised to construct and validate a prognostic model for osteosarcoma. The analysis involved the use of the LASSO Cox-regression method with the Glmnet R package to identify key PANoptosis-related genes. Differential gene expression analysis was conducted using the Limma R package, and model validation was performed using Kaplan-Meier survival analysis and time-dependent ROC curves.
Results: This model, derived from five key PANoptosis-related genes, demonstrated significant predictive capability for patient survival across training and validation cohorts. Further analysis confirmed the model's effectiveness and identified metastasis stage and risk scores as the robust independent prognostic indicators.
Conclusion: The prognostic model offers a novel tool for osteosarcoma prognostication and underscores the therapeutic potential of targeting PANoptosis-related pathways.
Key Words: PANoptosis-related genes, Osteosarcoma, Prognosis, Bioinformatics, Tumour micro-environment.
INTRODUCTION
Osteosarcoma is the predominant malignant tumour of bone, predominantly affecting adolescents and young adults. Globally, osteosarcoma accounts for approximately 3,500 to 4,000 new cases annually, with the highest incidence observed in adolescents and young adults.1 According to the Global Cancer Observatory (GLOBOCAN) 2020 data, the incidence rates vary geographically, with higher rates reported in regions such as Africa, Asia, and South America.2 Despite advancement in surgical techniques and adjuvant chemotherapies, the prognosis for patients with metastatic or recurrent osteosarcoma remains dismal.3,4 The five-year survival rate for these patients has stagnated at around 20% for the past several decades,5 highlighting an urgent need for the novel prognostic markers and therapeutic targets.
PANoptosis, a term recently coined to describe a form of programmed cell death that involves features of pyroptosis, apoptosis, and necroptosis,6 has emerged as a crucial mechanism in cancer biology. Contrary to classical pathways of cell demise, PANoptosis is orchestrated by unique triggers, molecular complexes, and executioners.7 The dysregulation of PANoptosis pathways has been implicated in the initiation, progression, and therapeutic resistance of various cancers, including colon and pancreatic cancer.8,9
Recent studies have shed light on the intricate roles of PANoptosis in the context of cancer, offering new insights into tumour biology and potential therapeutic strategies.10 Recent studies have demonstrated that the recruitment of Receptor-interacting serine / threonine-protein kinase 1 (RIPK1) to secondary signalling complexes, plays a critical role in the resistance of osteosarcoma to TRAIL-based therapies.11 Ning et al.’s work also revealed that osteosarcoma is characterised by the decreased expression of Caspase-8, which correlates with increased cell proliferation and reduced apoptosis, suggesting its roles as potential markers and therapeutic targets in osteosarcoma treatment strategies.12 Thus, understanding the role of PANoptosis-related genes in osteosarcoma could unveil new avenues for prognostication and treatment.
However, despite the potential of PANoptosis-related genes as prognostic markers and therapeutic targets in osteosarcoma, comprehensive models incorporating these genes are scarce. There is a critical gap in the literature regarding the systematic evaluation and validation of PANoptosis-related genes for predicting the prognosis of osteosarcoma patients. This study aimed to construct and validate a novel prognostic tool for osteosarcoma and highlight potential therapeutic targets based on PANoptosis-related genes using the LASSO-Cox regression.
METHODOLOGY
Transcriptomic data were extracted from the GEO database (GSE16088), which included samples from 6 normal individuals and 16 patients with osteosarcoma. This dataset served as the foundation for the investigation into the genetic alterations affecting PANoptosis-related genes within osteosarcoma.
For the training cohort, RNA-sequencing data alongside clinical details from 85 osteosarcoma patients were retrieved from the TARGET database.13 In contrast, the external validation cohort comprised RNA-sequencing data and corresponding clinical details for 53 osteosarcoma patients, sourced from the GEO database (GSE21257). Furthermore, the authors identified a set of 277 PANoptosis-related genes from an existing study,14 contributing to this study.
The GSE16088 dataset was leveraged to perform a comparative analysis on the expression levels of PANoptosis-related genes between osteosarcoma samples and normal tissue. The Limma R package was used to identify PANoptosis-related DEGs.15 The criteria for selection included genes exhibiting a |logFC| >0.585 and FDR <0.05. To identify differentially expressed PANoptosis-related genes with prognostic value, univariate Cox-regression analysis targeting overall survival (OS) was utilised.
To construct a prognostic model, TARGET-OS datasets were chosen as a training cohort, and PANoptosis-associated DEGs were integrated into a LASSO Cox-regression framework via the Glmnet R package.16 Subsequently, genes identified through LASSO regression were further analysed using multivariate Cox-regression with the survival package,17 aiming to develop a model capable of predicting patient outcomes. The risk score for each subject was computed based on the formula:
Risk score = Σ (expression (Genen) × coefficient (Genen) (expression (Genen) denotes the expression level of a specific gene, and coefficient (Genen) refers to its corresponding coefficient).
Based on the median risk score, patients were categorised into either low- or high-risk groups.
For external validation, the same prognostic model divided osteosarcoma patients from the GSE21257 dataset into comparable low- and high-risk categories, employing the identical median risk score criterion.
The independent predictive value of the signature was assessed using multivariate Cox-regression models. To evaluate the prognostic precision of the risk score across various datasets, Kaplan-Meier survival curves were utilised to assess the efficacy of the model in distinguishing between patient subtypes. The time-dependent receiver operating characteristic (ROC) curve and the incident / dynamic (I/D) area under the curve (AUC) metrics were employed, utilising the time ROC and Riskset ROC R packages, respectively.17,18
To ascertain the independent prognostic relevance of the risk score and assorted clinical parameters, both univariate and multivariate Cox-regression analyses were conducted. Factors exhibiting significance (p <0.05) in both analyses were identified as independent prognostic indicators and were subsequently incorporated into the development of a nomogram. The construction of this nomogram, aimed at predicting overall survival (OS), utilised the survival and rms R packages. The nomogram's predictive accuracy was quantified using the Concordance index (C-index) which varies between 0.5 and 1.0. A C-index approaching 1.0 denotes higher precision in the nomogram’s predictive ability. The nomogram’s performance was further assessed through Calibration curves for 1-year, 3-year, and 5-year OS predictions, comparing predicted outcomes with observed results. Ideally, a Calibration curve aligns closely with the 45-degree line, indicating accurate prediction.
The tumour immune microenvironment (TIME) characteristics across various molecular subtypes or PANscore categories were examined using the CIBERSORT algorithm. This involved assessing the infiltration of 22 immune cell types across all patient samples. Samples yielding p <0.05 were deemed accurately evaluated for immune cell infiltration and were selected for further analysis.
Furthermore, single nucleotide variation (SNV) and copy number variation (CNV) data from the TARGET-OS datasets were accessed. The mutational landscape was depicted using the waterfall function in the maftools R package, and Gistic 2.0 was utilised for CNV analysis.19
RESULTS
A comparative analysis of the PANoptosis-related gene expression between osteosarcoma samples and normal tissues, using the GSE16088 dataset, revealed significant differential expression. Employing the Limma R package, numerous PANoptosis-related differentially expressed genes (DEGs) meeting a stringent criteria of |logFC| >0.585 and FDR <0.05 were identified. These genes showed profound alterations in the osteosarcoma transcriptome, highlighting the potential mechanistic involvement of PANoptosis in tumourigenesis (Figure 1A), while their correlation coefficients were represented in Figure 1B.
Utilising the TARGET-OS dataset, a univariate Cox regression analysis of 277 PANoptosis-related genes identified 9 genes significantly associated with the survival of osteosarcoma patients (Figure 1C).
Figure 1: Analysis of differentially expressed PANoptosis-related genes and development of the prognostic model. (A) Heatmap of PANoptosis-related genes in the GSE16088 cohorts. (B) The correlation coefficient among the PANoptosis-related genes. (C) Forest plot of 9 PANoptosis-related genes significantly correlated with prognostic factors. (D, E) LASSO Cox-regression analysis was conducted on nine candidate PANoptosis-related genes to evaluate their prognostic significance, with the determination of the optimal penalty parameter for LASSO regression.To enhance the precision and simplify the predictive signature, LASSO regression and cox-stepwise regression analyses were employed to develop a prognostic model (Figure 1D and E), narrowing down the predictors to five essential genes: CRYBA1, CASP5, PIK3R5, CASP6, and PRKACB. The final prognostic model was established as follows:
Risk score = CRYBA1 × 0.0950938250121901 + CASP5 × 0.0613896146052882 - PIK3R5 × 0.51181272058199 - CASP6 × 0.103080159347395 - PRKACB × 0.2618011 75025256.
The prognostic model’s effectiveness was validated through Kaplan-Meier survival curve analysis in both the TARGET-OS and GSE21257 datasets. In the training set, the high-risk group as determined by the median risk score demonstrated significantly poorer survival than the low-risk group (p <0.001, Figure 2A and 2B). This distinction was consistently observed in the validation set, reinforcing the model’s predictive capability across the independent cohorts (Figure 2D and E) as well as the pooled dataset of TARGET-OS and GSE21257 (Figure 2G and H). Further analysis revealed that the prognostic risk score model had strong predictive capability for 1-year, 3-year, and 5-year patient survival rates, demonstrating excellent specificity and sensitivity in predicting survival outcomes (Figure 2C and F, and 3).
To evaluate the independent prognostic ability of the PANoptosis-related genes risk signature, both univariate and multivariate Cox-regression analyses were conducted, incorporating risk scores and various patient clinical characteristics. A univariate Cox-regression analysis (Figure 3A) identified metastasis stage and risk scores as significant independent prognostic factors (p <0.001). Subsequent multivariate analysis (Figure 3B) further delineated that among these, metastasis stage and risk scores consistently emerged as robust independent prognostic indicators (p <0.001). Leveraging the TARGET-OS dataset, a predictive nomogram (Figure 3C) integrating risk scores with clinical parameters was crafted, aiming to enhance the prognostication for patients. Combined with the calibration plots (Figure 3D), it was observed that the predicted OS was consistent with the actual observations. The nomogram’s efficacy was underscored through the ROC curve analysis (Figure 3E), revealing its superior predictive performance relative to the traditional clinical indicators. This superiority was further evidenced by the ROC curves for three- and five-year survival predictions (Figure 3E), showcasing the nomogram’s heightened predictive capability for patient outcomes.
This study explored the association between the PANoptosis-related genes risk score, tumour mutational landscapes, and immune infiltration profiles. The analysis revealed a distinct somatic mutation landscape across PANoptosis-related genes (PAN_DEGs) subgroups, highlighting the top 20 mutated genes including TP53, MUC16, and ATRX, among others (Figure 3F). Notably, TP53 emerged as the most frequently mutated gene, with a higher mutation rate observed in the low-risk score subgroup compared to the high-risk score group. Further, an immunological landscape analysis, conducted within both TARGET-OS (Figure 3G) and GSE21257 (Figure 3H) datasets, demonstrated significant variations in the distribution of 22 immune cell types across the subgroups. This comprehensive analysis underscores the intricate interplay between the PANoptosis-related genes risk score, tumour mutation burden, and immune cell infiltration, offering novel insights into the potential mechanisms driving these associations.
DISCUSSION
The findings of this study underscore the critical role of PANoptosis-related genes in the prognosis of osteosarcoma, highlighting a novel prognostic model that effectively predicts patient outcomes.
Osteosarcoma is linked to the dysregulation of cell death mechanisms and inflammatory responses. Genes associated with PANoptosis have been cited in the context of various other tumours.20 Members of the Caspase family are implicated in tumour progression.21 The identified genes within the model may serve as potential targets for the novel therapeutic strategies, aiming to modulate the pathways of programmed cell death to improve patient outcomes.
This study contributes to this growing body of knowledge by exploring the correlation between tumour immune infiltration levels, genetic changes, and risk stratification in cancer patients.22
Research has demonstrated that the tumour microenvironment, characterised by the immune cell infiltration, plays a pivotal role in cancer progression and therapeutic response.23 Tumour mutational burden (TMB) significantly correlates with immune scores and tumour purity, indicating that high TMB levels are associated with altered immune infiltration​. Differences in the proportion of various types of immune cells, such as CD8+ T cells and macrophages, were linked to survival outcomes.24 A study on prostate cancer reinforced the connection between genetic variations, immune cell infiltration, and patient risk stratification. It was found that certain somatic mutations were more prevalent in the high-risk groups, which also correlated with differences in tumour-infiltrating immune cells.25
This research aligns with these findings, suggesting that the interplay between genetic variations and the level of tumour immune infiltration can significantly impact the prognosis of cancer patients. This underscores the need for further investigation into the mechanisms driving these associations and their implications for cancer treatment and patient care.
To enhance the generalisability and clinical applicability of this prognostic model, future studies should focus on validating the model across larger, more diverse populations, and in various clinical settings. The lack of prospective validation is a limitation of this study, as real-time patient data would provide a more accurate assessment of the model’s performance. The authors recommend that future research should include prospective studies to confirm the utility of PANoptosis-related gene signatures in osteosarcoma prognosis and to explore their potential in guiding therapeutic strategies.
The Kaplan-Meier curves demonstrate a clear separation between the high-risk and low-risk patient groups, indicating the model’s ability to stratify the patients based on survival probability.
Furthermore, the time-dependent ROC curves underline the predictive precision of the risk score, with AUC values exceeding 0.70 across different time points (1-year, 3-year, and 5-year survival), reinforcing the model’s robustness. These graphical insights confirm that the proposed model not only offers a novel tool for osteosarcoma prognostication but also holds potential for guiding therapeutic decisions by identifying high-risk patients who may benefit from more aggressive treatment strategies.
The integration of PANoptosis-related gene analysis in osteosarcoma prognostication, supported by rigorous graphical validation, highlights the clinical utility of this model. Future research should focus on prospective validation in diverse populations to further solidify the model’s applicability and explore its potential in guiding personalised treatment approaches.
CONCLUSION
This study elucidates the prognostic significance of PANoptosis-related genes in osteosarcoma by presenting a validated model that effectively predicts patient outcomes. The graphical representations, including Kaplan-meier survival curves and ROC curves, provide strong visual evidence of the model’s discriminatory power.
ETHICAL APPROVAL:
The study was initiated after obtaining ethical approval from the Ethical Review Committee of the Second Xiangya Hospital, Hunan, Xiangya, China (Approval No.: LYEC2024- 0093).
PATIENTS’ CONSENT:
This study did not involve direct patient data; hence, no informed consent was required.
COMPETING INTEREST:
The authors declared that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
AUTHORS’ CONTRIBUTION:
XP: Conducted the investigation and drafted the manu-script.
XP, WW: Conceived and designed the study, critically revised the manuscript for the important intellectual content.
Both authors approved the final version of the manuscript to be published.
REFERENCES